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Abstract: We report on our computation of the leading hadronic contribution to the anoma- 
lous magnetic moment of the muon using two dynamical flavours of non-perturbatively 
0(a) improved Wilson fermions. The strange quark is introduced in the quenched approxima- 
tion. Partially twisted boundary conditions are applied to improve the momentum resolution 
in the relevant integral. Our results, obtained at three different values of the lattice spacing, 
allow for a preliminary study of discretization effects. We explore a wide range of lattice 
volumes, namely 2 fm < L < 3 fm, with pion masses from 600 to 280 MeV and discuss 
different chiral extrapolations to the physical point. We observe a non-trivial dependence of 
a^^^ on ttZtt especially for small pion masses. The final result, a^^^ = 618(64) x 10"-'^'^, is 
obtained by considering only the quark connected contribution to the vacuum polarization. 
We present a detailed analysis of systematic errors and discuss how they can be reduced in 
future simulations. 
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1 INTRODUCTION 



The magnetic moment of a charged lepton is extracted from the vertex function describing the 
interaction between the lepton and a photon in the limit of vanishing photon momentum. The 
corresponding anomalous magnetic moment ai is then defined as half the difference between 
the gyromagnetic factor g and its classical value of 2, i. e. ai = [gi — 2)/2. In the case of the 
electron the quantity is dominated by QED contributions. The one-loop result was obtained 
by Schwinger more than sixty years ago [1] , and since then Ce has reached an accuracy better 
than one part per billion on both the theoretical and the experimental sides, which yield 
results in beautiful agreement (see section 3 of [2], and references therein). 

The anomalous magnetic moment mediates helicity flip transitions [3], which implies 
that quantum corrections due to heavier particles of mass M, in the Standard Model or 
beyond, are proportional to m^/M^. For this reason the muon anomalous magnetic moment 
is regarded as a sensitive probe for effects of nearby New Physics. However, by the same 
argument, given that m^j < m-Tr, the hadronic contributions to are large and notoriously 
difficult to quantify. While the experimental and theoretical estimates have each reached 
similar levels of precision of 0.5 ppm, a tension by 2 or 3 standard deviations between theory 
and experiment persists [4, 5]. Before invoking New Physics as the reason for this tension 
the theoretical result and, in particular, all contributions due to hadronic effects, must be 



phenomenological approach based on the evaluation of a dispersion integral. In the low-energy 
regime the spectral function in the integrand must be determined experimentally, either from 
the cross section e'^e~ — ?■ hadrons or from the rate of hadronic r-decays. Both methods 
suffer from different systematics [2, 5] and yield results in slight tension among each other. 
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None of them however reduces the discrepancy between theory and experiment on below 
the two standard deviation level. A purely theoretical estimate of aj^^*-* from a first-principles 
approach is clearly desirable, and the present work represents our first step in that direction 
by using the lattice regularization of QCD. 

The hadronic vacuum polarization contribution to has received considerable atten- 
tion from the lattice community in recent years. Initial studies have been performed in 
the quenched approximation [6, 7] and in the theory with two [8] and three dynamical 
fiavours [9, 10]. Compared to the determination of "standard" quantities such as hadron 
masses, quark masses and decay constants (see [11] for a review), lattice calculations of a^^'~' 
are extremely challenging. The relevant lattice quantity (the hadronic vacuum polarization 
discussed in the next section) receives contributions from quark disconnected diagrams, which 
are intrinsically noisy and difficult to estimate with good statistical accuracy at a reasonable 
numerical cost. Moreover, the dependence of the hadronic vacuum polarization on the mo- 
mentum transfer must be accurately traced down to momenta of order and beyond. This 
value is well below the lowest Fourier momentum (27r/L)^ which can be reached in current 
lattice QCD simulations. As a consequence finite size effects may conceivably be large on 
results obtained within the "standard" approach. In addition, the vacuum polarization re- 
ceives sizeable contributions from the low-lying vector resonances. Those should be properly 
accounted for in simulations at sufficiently light quark masses for the lowest vector meson to 
be a resonance and including the full dynamics of the strange and charm quarks. Finally, as 
the contribution to a^^^ from isospin breaking effects may be of the same size of its present 
uncertainty, those will have to be included in lattice computations (possibly along the lines 
discussed in [12-14]) for them to have a crucial impact on (g — 2)^ phenomenology. 

In [15, 16] we have shown how (partially) twisted boundary conditions [17-19] can be 
used to improve the momentum resolution in the connected part of the hadronic vacuum 
polarization, and we have obtained an estimate of the disconnected contribution in Chiral 
Perturbation Theory, thus addressing the first two systematic effects discussed above. Here we 
report on numerical results obtained in that setup and show how the fits to the momentum de- 
pendence of the vacuum polarization function get substantially stabilized in this way, thereby 
improving the accuracy of our estimates. Partial results have already appeared in [20, 21]. 
In addition we present a thorough investigation of all sources of systematic errors, including 
the modelling of the g^-dependence of the vacuum polarization, chiral extrapolations, lattice 
artifacts and finite volume effects. 



2 DEFINITIONS AND LATTICE SETUP 



The Euclidean hadronic vacuum polarization (VP) tensor is defined as 
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where (x) = ^ Q ftp jr{x)j^'ipf{x). For N{ = 2 the quark fields ipf{x) are taken from the 
/=i 

set (ipujipd) with Qf = (2/3,-1/3) denoting the electric charges in units of the elementary 
one. For A'^f = 2 + 1, / = {u,d,s) and Qf = (2/3,-1/3,-1/3). Euclidean invariance and 
current conservation imply 

^lf\Q) = i9,uq' - q,qu)n(''^\q') • (2.2) 

For space-like momenta, the relation between I[^f^^\q'^) and the lowest order hadronic contri- 
bution a^^'^ to the anomalous magnetic moment of the muon has been derived in [6, 7, 22] 
and reads (suppressing the index A^f) 

,HLO ^'^ 



.IT 

where 

m 



«r = F / dq^fiq^Mq^), (2.3) 



Wz^n - q^z) q^-Jq^ + ^m^q^ 

1 + mj^q'^Z'^ Im'jj^q'^ 



and n(g2) =47r2 [n(g2) _ n(0)] . 
2.1 Lattice regularization 

We perform our computation on a subset of the gauge configurations generated within the CLS 
initiative [23] for two flavours of non-perturbatively 0(a) improved Wilson fermions [24] and 
using the standard Wilson plaquette gauge action. The simulation parameters are collected 
in Table 1, including the values of the twist angle 6, whose role will be explained below. 
Following [7], we have implemented the one-point-split conserved vector current 

V^{x) = \ {i^f[x + a/i)(l + 7^)C^^(x) i^f[x) - V^^(x)(l - i^)U^{x) 4^f{x + afi)) , (2.5) 

where £ SU{3) represents the gauge link in the positive fi direction, and / is again a 
flavour index. The lattice version of the vacuum polarization tensor then reads 

Nf Ni 

KfHx) = a\j^{QjVfix)) J^iQf'V/m) . (2.6) 
/=i /'=i 

By Fourier transforming the expression above into momentum space one gets 

n(^*)(g) = E e^'?(-+"'^/2--/2)n2:')(^) , (2.7) 

where q^ = 27rn^/L and q^ = ^ sin (^) with G 0, 1, . . . L/a — 1. We restrict our attention 
to the case ^ v to avoid mixings of the composite field V^(2;)VJ/(0) with lower dimensional 
ones in the limit x — t- 0. It then follows that Il^if^^\q) fulfils the Ward identities 



4^')(^)=nW)(g)g. = 0. (2.8) 
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^ meas. 


A3 


64 


32^ 


5.20 


471 


0.079 


0.8 


1.8 


2.6 




532 


A4 


64 


32^ 


5.20 


362 


0.079 


0.8 


1.8 


2.6 




800 


A5 


64 


32^ 


5.20 


317 


0.079 


0.8 


1.8 


2.6 




432 


E4 


64 


32^ 


5.30 


601 


0.063 


0.8 


1.8 


2.6 


0.13605 


648 


E5 


64 


32^ 


5.30 


447 


0.063 


0.8 


1.8 


2.6 


0.13574 


672 


F6 


96 


48^ 


5.30 


324 


0.063 


0.4 


1.9 


2.3 


0.13575 


804 


F7 


96 


48^ 


5.30 


277 


0.063 


0.4 


1.9 


2.3 


0.13570 


820 


N4 


96 


48^ 


5.50 


541 


0.050 


0.8 


1.9 


2.6 


0.13639 


532 


N5 


96 


48^ 


5.50 


431 


0.050 


0.8 


1.9 


2.6 


0.13629 


644 



Table 1. Summary of simulation parameters. Measurements are performed on configurations sepa- 
rated by 8 units of molecular dynamics time at least, where one molecular dynamic unit corresponds 
to an evolution by r = 1 in Hybrid Monte Carlo algorithms. The lattice spacings and the values of 
the hopping parameter Kg corresponding to the strange quark mass are taken from [25] and have been 
determined using the definitions and some of the results in [26-28] . The values of the lattice scale must 
be regarded as preliminary. The masses of the lightest pseudoscalar state are taken from [29] . The 
twist angle 6 is applied in the spatial x-direction only. 



Consequently on the lattice we can relate the scalar vacuum polarization to the tensor one, 
by mimicking the continuum relation 

^if\Q) = iQ\u - q,qu)n^'''H f) • (2.9) 

The scalar vacuum polarization extracted in this way approaches its continuum counterpart 
with a rate proportional to a, where the 0(a) effects appear due to off-shell contributions in 

After performing the Wick contractions in Illf^^\x) one realizes that connected as well 
as disconnected quark diagrams contribute. We neglect here the disconnected terms, as 
done in most of the existing lattice computations in the literature. In [8] such contributions 
were included for almost half of the ensembles used and found to be negligible within the 
quoted errors. Following the arguments in [16], the disconnected diagrams are expected to 
decrease the value of f[^f^^\q) by about 10%, which demonstrates the limited accuracy of 
lattice calculations without their proper inclusion. 

2.2 Twisting the connected part of the vacuum polarization 

By modifying the spatial boundary conditions on the quark fields entering the vector current it 
is possible to improve the momentum resolution in Il^'^^\q'^) and to access momenta different 
from the integer multiples of 2tt/L. Namely, imposing the condition 

^{x + Lk) = e''^>=ip{x) (2.10) 
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is equivalent to boosting the momenta in the quark propagator by 9k/L. This technique 
can be used to modify and refine the lattice dispersion relation of, for example, charged 
pseudoscalar mesons, or of any momentum dependent quantity, as long as there are no strong 
final state interactions [18]. This remains true, up to exponentially suppressed finite size 
effects, for partially twisted boundary conditions [17-19], where only (some of) the valence 
quarks satisfy twisted boundary conditions, while the sea quarks fulfil periodic boundary 
conditions. Note, however, that the effect of partial twisting obviously vanishes whenever it 
is applied to flavour-singlet quantities, as the vector current discussed above. 

We have shown in [16] that, by introducing a sufficiently large number of valence quarks 
degenerate with the u, d . . . flavours, the quark connected and disconnected parts of the cor- 
relator Il\J^^\x) can be rewritten as independent correlation functions in Partially Quenched 
QCD. In such unphysical theories each quark diagram has an unambiguous field-theoretic 
expression and well-defined continuum and infinite volume limits. The formulation naturally 
turns the connected contribution into a correlator of flavour-off-diagonal vector currents and 
thus twisting can be applied to induce arbitrary momentum. This amounts to a simple mod- 
ification of Eq. 2.7 (restricted to its connected part) in which ^ qk — Ok/ L and q changes 
accordingly. 

In practice we have applied twisting in the x-direction only and remained with periodic 
boundary conditions in the other directions. The twist angles in Table 1 have been chosen in 
order to achieve a fixed and large density of q^ values between consecutive Fourier modes and 
in order to reach a lowest non-zero around m^. For the Fourier modes we have considered 
all the integers values between (0, 0, 0, 0) and (2, 2, 2, 2) in units of 27r/L. 



we introduce a maximum momentum ^^^x adopt several different fit ansdtze for the low 
q"^ region below Q^ax- We additionally impose a matching condition with perturbation theory 
at (f = q^i^^ and hence use the perturbative expression to describe the high q'^ region. 

The results for a^^'~' at different values of the mass of the lowest pseudoscalar state 
are then extrapolated to the physical point using a functional form inspired by Chiral 
Perturbation Theory. 

3.1 Fitting procedure 

In order to check for the stability of our results against variations in the fitting procedure we 
use three different ansdtze to describe H^^^^q"^) in the region 

A) a model-independent (2, 3) Pade ansatz defined by a degree 2 over a degree 3 polynomial 
in f (ricoeff = 6), i.e. 



3 RESULTS 



In order to compute we must obtain a continuous description of Il^^^\q'^) to determine 
n(^4)(0) and fl^^^^q"^) entering the integral in Eq. 2.3. At each value of the lattice spacing a 
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B) a vector-dominance model with a single pole (ncoefr = 3), i.e. 

n(N^\f)=a+^^, (3.2) 
+ 

C) a vector dominance-model with two poles^. This form {ricoeS = 5) has been suggested 
and used in [10] and reads 

n(^f)(^2) = a + + , (3.3) 

subject to the additional constraint < e^. 

We perform correlated fits. For these to be reliable and to produce a meaningful correlated 
large statistics is needed. In [30] a toy model was used to provide the thumb-rule D'^, 
where N is the number of measurements and D the number of degrees of freedom in the fit. 
We therefore randomly select about 25 values of (f and indeed observe stable correlated fits 
as well as the absence of numerical problems when inverting the average covariance matrix. 
The errors on the fit parameters are estimated and propagated to a^^'~' by performing the 
fits and the integral for each bootstrap sample. For the fit ansatz C we compare in Fig. 1 the 
value of the smaller mass parameter obtained from the fit with the mass measured from the 
exponential decay of the vector two-point function. We call the latter naive vector mass, as 
it assumes the lightest vector meson to be stable at all values of m.j^ considered. While this 
may not be the case for our most chiral points, for the sake of this qualitative comparison 
such a definition seems nonetheless adequate. A separate study will have to be devoted to 
the measurements of the vector resonance parameters following [31, 32]. The figure shows 
overall reasonable agreement, with discrepancies which stretch to 15% at most. 

As a constraint on the high-^^ region we always include perturbation theory in the MS 
scheme at the scale 2 GeV, using the leading-order expression in as for n(^f)(g2) from [33]. 
Although the NLO in is available from the same reference, for our purposes it is enough to 
include the leading-order, where in addition no internal fermionic loops appear and therefore it 
is rather simple to apply the formulae to the case of two light dynamical quarks and a quenched 
strange. In addition, as it will become clear in the following, the perturbative contribution 
to a^^*-* estimated in this way is extremely small. To evaluate the perturbative formula 
we use the non-perturbative Aj^ parameter for = 2 from [34] and the non-perturbative 
renormalization factors in [35-37] to relate the lattice quark masses to their values in the 
MS scheme at 2 GeV. As discussed in [7] the function computed on the lattice and 

in the continuum dimensional regularization can differ by an un-physical and non-universal 
integration constant. We fix the value of this scheme dependent constant by matching the 
perturbative expression for Il{q'^) to our non-perturbative data at = g^ax- Moreover the 
perturbative prediction diverges as — )• and therefore even after this step perturbation 

'^Up to a re-definition of the fitting parameters the ansdtze B and C are equivalent to Fade approximations 
of degree (1,1) and (2,2) respectively. 
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Figure 1. Comparison between the smaller mass parameter from fit ansatz C and the vector mass 
measured from the vector two-point function as a function of m^. 



theory can be expected to provide a good description of our data only at sufficiently large 
values of g^. We require the resulting function not only to be smooth but also to have a 
smooth first derivative on the whole real axis, and in particular at cf' = g^ax- Hence, the 
matching produces one non-trivial relation among the parameters in the fit ansdtze above, 
reducing their number ripar from npar = '^coeflf to ncocff — 1- 

The fit ansatz B always produces reduced larger than 5 and we therefore discard it. 
We monitor the stability of the fit results for a|^^° in the way shown in the left panel of Fig. 2 
where, as an example, the values obtained from fits of type A and C to the measurements 
performed on the F6 ensemble are plotted. A cut x^/dof < 2.5 is applied in producing the 
figure. The two functions give compatible results for ^max~ 1-5 GeV^ and in the end we opt 
for the results from ansatz A on all the ensembles, as this yields in most cases the smallest 
values and the smallest errors on aj^^^ . The values of (/^ax chosen differently for different 
ensembles as the g^-ranges explored depend on the lattice spacing and on L/a. Specifically 
the maximum momentum lies in the range 2.4 GeV^ ^ ^max ~ 5.5 GeV^. In order to quantify 
the improvement brought by the use of partially twisted boundary conditions we repeat in 
Fig. 2, right panel, the comparison between fit ansdtze A and C on the F6 ensemble but 
restricting the data to the Fourier modes only. Twisting clearly helps in stabilizing the fits 
and the plateau values are reached "earlier" in ^^ax- This increases our confidence in the 
fitting adopted procedure. Statistical errors are also slightly smaller in the left panel of Fig. 2 
with respect to those in the right one. 
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Figure 2. Left: Fit results from the F6 ensemble using the ansdtze A and C in the text, plotted as 
a function of g,nax- R-ight: Equivalent plot obtained by using the Fourier momenta only. 



We summarize our results on for the two-flavour theory and the case with an 

additional quenched strange quark in Tables 2 and 3 respectively. 
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274.2(9.8) 


0.96 


282.4(9.1) 


0.97 


A4 
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Table 2. Results for a^^o in the iVf 



2 theory. 



How statistical errors of order 5-10% come about can be easily understood from Fig. 3. 
In the left panel we show the result for n*^^^(g^) together with the fit function (ansatz A) 
and the perturbative curve for the F6 ensemble. The horizontal axis is divided into four 
regions. Region I, is defined by < < (not visible on the left panel) where the q^- 
dependence is only constrained by smoothness requirements on the fit function. There is 
no direct measurement of the vacuum polarization function here, although we could have 
tuned the 0-angle in order to penetrate this region. However, the error on the corresponding 
value of the integrand in Eq. 2.3 would have turned out to be larger than 100%. Region 
II, where < < 0.2 GeV^, is accessible only thanks to the use of partially twisted 
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Table 3. Results for a^^° in the theory with two dynamical light quarks and a quenched strange 
quark. 



boundary conditions; 0.2 GeV^ indeed coincides witfi tlie smallest Fourier mode on the F6 
ensemble. Region III, 0.2 GeV^ <q^< gLx, is the region which is accessible here by Fourier 
modes. Region IV, > ^^ax; is the "perturbative" region. The right panel shows the 
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Figure 3. Left: The subtracted vacuum polarization &^\q-^) computed on the F6 ensemble. The 
blue solid is the result from fit ansatz A matched to two-loop perturbation theory (light blue line). 
Right: The different contributions to a^^'~' broken down to different momentum ranges (see text), 
indicated in both panels by corresponding colours. 



relative contributions from the different regions to the integral in Eq. 2.3. Since region I 
contributes about 25% and is not constrained by any direct measurement of n^^)(^^), its 
uncertainty dominates the overall statistical error on a|J^^, with an ambiguity which, within 
this approach, inevitably amounts to about 5%. The uncertainty on the contribution from 
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region II is of the same order but slightly smaller, whereas the uncertainties on regions III 
and IV are sub-dominant. 

3.2 Chiral extrapolations 

We show our results for the N{ = 2 and the Nf = 2 plus a quenched strange quark 

theories as functions of in Fig. 4. The blue curves represent our preferred chiral fits and 
we will discuss in the following how those have been obtained. We start from the observation 
that in the mass-region explored a curvature is clearly visible in our data. Indeed a linear 
fit in that region would fail in producing acceptable x^/dof values. We take the observed 
curvature into account in three different ways. 

• We perform a chiral fit inspired by the functional forms derived at NLO in for 
different well known chiral corrections. These include a chiral logarithm, which should 
account for the curvature. 

• We re-analyze our data following the alternative procedure presented in [8], which is 
expected to remove the curvature due to the lowest lying vector resonance. 

• We impose a cut on the pion mass (< 400 MeV) and linearly extrapolate our data at 
the four lightest values of m^r in the Nf = 2 theory. In doing so we negelect cutoff effects 
and fit the /3 = 5.2 and /? = 5.3 results simultaneously. Such a procedure produces a 
value for a^^^ at the physical point which is well consistent with the result from the 
xPT-inspired fit ansatz in Eq. 3.5 although slightly lower (namely a^^^ = 508 (62)). 

We do not further discuss the third analysis and present in detail the approaches used for the 
first two fitting procedures. 

Since little is known about the functional form describing the dependence of a^^^ on 
the light quark masses, we conservatively adopt a form inspired by xPT to extrapolate our 
data to the physical point. Namely we use the fit function 

aHLO(^^) = aHLO(o) + B {am^f + C (am^f Hiam^f) , (3.4) 

with a^'^^, B and C as free parameters. Such a functional form can be derived for the 
connected part of a^^'^ from the expressions for the (connected) vacuum-polarization obtained 
in [16] at NLO in xPT. In addition, since our current set of ensembles covers a wide range 
of pion masses at only one value of the lattice spacing (a = 0.063 fm), and in order to avoid 
mixing of cutoff effects with chiral effects, we decide to extrapolate the /3 = 5.3 data only. The 
other data sets in Table 1 are used to asses lattice artifacts and finite size effects. Eventually, 
extrapolations to the continuum limit at fixed pion masses will have to be performed before 
the chiral extrapolation. 

The resulting curves are shown in Fig. 4. We find that the function describes the whole 
set of data points quite well, even those which are not included in the fit. In the two-flavour 
case we obtain at the physical point 

a|^LO ^ (g2) . iQ-io [Ni = 2] , (3.5) 
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Figure 4. Left: The simulation results for the hadronic contribution to computed using two 
flavours, shown as a function of m^. The chiral extrapolation (blue curve) is performed for the 
ensembles at /? = 5.3 using an ansatz motivated by chiral perturbation theory. Right: Corresponding 
results for a}^^^ including a quenched strange quark. The value from the PDG [38] is also shown for 
illustration. 



which is well consistent with the recent Ni = 2 result in [8], whereas the inclusion of a 
quenched strange quark gives 

a|JLO ^ 618 (58) • 10"^° [iVf = 2 + Ij , (3.6) 

where the errors are purely statistical. 

In [8] an alternative extrapolation procedure was proposed, with the aim of reducing 
the TTiTT-dependence. It can be motivated starting from a vector dominance description of 
n(^f)(^2), i.e. 



m\e,m^) oc g\ 2 , \ ^ -2 ' (3-7) 

where the m7r-dependence is explicitly indicated. The quantity gy is related to the vector 
decay constant /v by gy = fy/my. The dependence of gy on is neglected based on the 
numerical observation in Fig. 1 of [39]. Eq. 3.7 then leads to 



2 

mr, 



dq^f{q')m\q^m^) « . (3.8) 



It is easy to see that, by rescaling the argument of the function / from to hq'^, one obtains 

Under these assumptions the choice h = ^ in [8] would therefore remove the dependence 

on m-jr (up to the mild one in ^v) producing in addition the right physical result since /i — )■ 1 
in that limit. Such a rescaling has the further advantage of making the dimensionless quantity 
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a^^f* independent from the lattice spacing [8], as opposite to a^^^ where the lattice spacing 
value is needed to convert in the weight-function / (see Eq. 2.4) to lattice units. 

In order to estimate the uncertainty in our chiral extrapolation we perform the analysis 
described above, using the naive vector mass to define h. Results for are shown in Fig. 5. 
We expect from the observation in the previous section on the high x^/dof values produced by 



Nf = 2 



Nf = 2+l„ 




P =5!2,L = 2.4iiu I 
P =5.3,L = 2.0fm 
5.3,L = 3.0fm 
P =5.5,L = 2.5fm 
linear fit 



0.05 0.1 0.1-5 0.2 0.25 0.3 0.35 0.4 




P = 5.3. L = 2.0 fm 
P = 5.3.L = 3.0fm 
P = 5.5. L = 2.5 fm 
PDG 
linear fit 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
IGeV'l 



Left: Nf = 2. Right: 



Figure 5. Results and chiral extrapolations for with h = 

iVf = 2 + Iq. In both cases the leftmost value shows the result from the corresponding extrapolation 
in Fig. 4. 

the fit ansatz B (single vector dominance) that some dependence on will remain. However, 
the rescaling clearly renders such a dependence mild, its effect being particularly strong for 
heavy pion masses. The curvature visible in Fig. 4 has almost completely disappeared and 
we therefore linearly extrapolate the results from the three smallest pion masses at f3 = 5.3 
to the physical point. There is no obvious theoretical reason why a non-linear dependence 
could not survive, however, we do not have a large sensitivity to those terms in the mir-range 
explored here. We obtain in this way the following results at the physical point 



,HLO 



550 (42) • 10"^° [Nf = 2] 



and 



HLO 



626(38) • 10 



-10 



[Ni = 2 + Ic 



(3.10) 



(3.11) 



which are consistent with those in Eqs. 3.5 and 3.6. Notice that, as a consequence of the 
rather moderate chiral extrapolation performed, they are also compatible with the values of 



,HLO 



directly estimated at our most chiral point (F7 ensemble). 



3.3 Residual cutoff and finite size effects 

The set of our simulations covers a rather wide range of lattice volumes and lattice spacings. 
We are therefore able to address the issue of finite volume effects and cutoff effects. This 
study is not yet complete and should be extended to lower pion masses. 
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3.3.1 Cutoff effects 

In general we expect 0(a) discretization effects in a^^^ as the vacuum polarization receives 
off-shell contributions and the lattice regularization used here is only on-shell improved. Some 
indications on the size of cutoff effects can be gathered from Fig. 4, where at least for low 
pion masses discretization effects appear to be rather small and below our statistical errors. 
However, we believe it is more instructive to look directly at lattice artifacts in Il^^^\q'^), as 
those conceivably depend on g^. In Fig. 6 we compare the subtracted vacuum polarization 
from the ensembles N5 and A3, which have rather similar spatial extensions and pion masses 
but different values of the lattice spacing. As expected, lattice artifacts mainly distort the 



fTl2)(q2) 




Figure 6. The subtracted vacuum polarization from two ensembles with roughly the same pion 
mass m^r — 450 MeV and size L ~ 2.5 fm, but different lattice resolutions {a{A3) — 0.079 fm vs. 
aim) = 0.05 fm). 

function at large values of cf', in the region which contributes little to aj^^*~* (see Fig. 3). 
3.3.2 Finite size effects 

As for the case of cutoff effects we look at finite size effects directly in fl^^^\(f'), since 
different states are expected to contribute to the vacuum polarization when cf' is varied (see 
also [40] ) . In Fig. 7 (left panel) we compare the results from the E5 ensemble with those from 
a simulation at the same bare parameters but different L/a (D5 ensemble, not in Table 1), 
specifically L/a = 24 instead of 32. Finite size effects are clearly visible in this rather extreme 
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setup (niTjL ~ 3.4 for D5), but they seem to drop below our present statistical error as L is 
made larger, as the right panel of Fig. 7 shows. Note, however, that in this second comparison 
both L and a differ. 




Figure 7. Comparison of results for tl{q^) at different values of L. Left: The lattice E5 {L/a = 32, 
a — 0.063 fm) is compared to a L/a = 24 lattice with the same resolution (lattice D5, not used here 
in the chiral extrapolation). Right: F6 {L/a = 48, a = 0.063 fm) vs. A5 {L/a = 32, a = 0.079 fm). 

4 CONCLUSIONS 

We have presented a calculation of the leading hadronic contribution to the anomalous mag- 
netic moment of the muon on the lattice with two dynamical flavours and a quenched strange 
quark. We have discussed technical improvements, which led to a better determination of 
the external momentum dependence of the vacuum polarization Il^^^\q'^). Specifically, to- 
gether with [41], this paper contains on of the first numerical applications of partially twisted 
boundary conditions to a quantity containing flavour-singlet currents. The approach follows 
from the theoretical setup devised in [16] . We restrict the computation to the connected part 
of &'^^\q'^) only. As it is clear that presently the accuracy of lattice results is not yet at the 
level required by phenomenology, the main goals of this paper are a precise assessment of the 
different sources of uncertainty and an estimate of their size. The main results are listed in 
Eqs. 3.5 and 3.6, where statistical errors only are included. In the following we estimate the 
uncertainties due to the modelling of the g^-dependence of the vacuum polarization, chiral 
extrapolations, lattice artifacts and finite volume effects. These uncertainties are quantified 
by changing one ingredient at a time in the conservative analysis procedure, which we have 
followed in producing our main results. 

In detail, on top of the statistical error we identify the following sources: 

• Fitting procedure. We repeat the entire analysis by using fit ansatz C instead of 
ansatz A everywhere and obtain 

aJ^LO = (55) . iq-w ^ 3] , [ansatz C] , (4.1) 
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a^^^ = 615 (56) • 10"^" [iVf = 2 + Iq] , [ansatz C] , (4.2) 
which suggests an uncertainty weU below our statistical errors. 

• Chiral extrapolation. By comparing the results in Eqs. 3.5 and 3.6 to those in 
Eqs. 3.10 and 3.11, we conclude that at the moment this systematic cannot be resolved 
with our statistical errors. 

• Cutoff and finite size effects. As discussed in Sect. 3.3 these effects appear to be 
small at the volumes, masses and lattice spacings considered here, but a more compre- 
hensive study is required. 

• Uncertainty in the lattice spacing. In [25] the lattice spacing at /3 = 5.3 is given 
as a = 0.063(3) fm, by combining statistical and systematic errors in quadrature. By 
repeating the analysis of our (/3 = 5.3) data using a = 0.066 fm instead of a = 0.063 fm 
we arrive at 

flHLO ^ ^Qg) . ^Q-io [iVf = 2] , [a = 0.066 fm] , (4.3) 
^HLO ^ Qj-^ (g4) . ^Q-io [^Tj = 2 + Iq] , [a = 0.066 fm] . (4.4) 

We will include this systematic by taking half the difference with respect to the results 
in Eqs. 3.5 and 3.6. 

In general we see that at the moment most of these systematic effects, with the exception of 
the uncertainty on a, are well below our statistical errors. They will become relevant and will 
have to be more precisely estimated once the latter are reduced. 

We quote as final results the values determined at f3 = 5.3, having used the other two 
lattice spacings to estimate systematic uncertainties. No continuum extrapolation has been 
performed yet. The present computation of the connected contribution to the anomalous 
magnetic moment of the muon then gives 

a^LO ^ (gg) . iQ-io ^ 2] , (4.5) 

and 

a|^LO ^ g^g (g4) . ^Q-io [Ni = 2 + Iq] , (4.6) 

where we have combined in quadrature statistical and systematic errors. 

The overall error can definitely be reduced, as it is still statistics-dominated. However, 
once it is evaluated, the contribution from disconnected diagrams, which is estimated to be 
around 10% (see [16]), will become the main uncertainty. We therefore consider studying the 
accuracy that can be reached on the numerical estimates of the disconnected contribution as 
a priority. A combination of the methods discussed in [42-45] seems very promising in this 
respect, and we have also started implementing similar techniques in the context of mesonic 
three-point functions [46] . Further improvements include considering lighter pion masses and 
enlarging the set of simulations at /3 = 5.5. 



15 



Acknowledgments: We thank Fred Jegerlehner, Rainer Sommer, Gilberto Colangelo, Achim 
Denig, and Harvey B. Meyer for useful discussions. We are grateful to our colleagues within 
the CLS project for sharing gauge ensembles. Calculations of correlation functions were 
performed on the dedicated QCD platform "Wilson" at the Institute for Nuclear Physics, 
University of Mainz. We thank Dalibor Djukanovic for technical support. This work was 
supported by DFG (SFB443) and the Research Center EMG funded by Forschungsinitiative 
Rheinland-Pfalz. 



References 

J. S. Schwinger, On Quantum electrodynamics and the magnetic moment of the electron, Phys. 
Rev. 73 (1948) 416-417. 

F. Jegerlehner and A. Nyffeler, The Muon g-2, Phys. Kept. 477 (2009) 1-110, 
[arXiv: 0902. 3360]. 

W. S. Cowland, On Schwinger's theory of the muon, Nucl. Phys. 8 (1958) 397-401. 

M. Davier, A. Hoecker, B. Malaescu, and Z. Zhang, Reevaluation of the Hadronic Contributions 
to the Muon g-2 and to a{M'^), Eur.Phys.J. C71 (2011) 1515, [arXiv : 1010 .4180]. 

F. Jegerlehner and R. Szafron, p*^ — 7 mixing in the neutral channel pion form factor |i^7rP 0''>T'd 
its role in comparing e'^e~with r spectral functions, Eur. Phys. J. C71 (2011) 1632, 
[arXiv: 1101.2872]. 

T. Blum, Lattice calculation of the lowest order hadronic contribution to the muon anomalous 
magnetic moment, Phys. Rev. Lett. 91 (2003) 052001, [hep-lat/0212018]. 

QCDSF Collaboration Collaboration, M. Gockeler et. al.. Vacuum polarization and hadronic 
contribution to muon g-2 from lattice QCD, Nucl. Phys. B688 (2004) 135-164, 
[hep-lat/0312032]. 

X. Feng, K. Jansen, M. Petschlies, and D. B. Renner, Two-flavor QCD correction to lepton 
magnetic moments at leading-order in the electromagnetic coupling, arXiv: 1103.4818. 

C. Aubin and T. Blum, Calculating the hadronic vacuum polarization and leading hadronic 
contribution to the muon anomalous magnetic moment with improved staggered quarks, 
Phys. Rev. D75 (2007) 114502, [hep-lat/0608011]. 

P. Boyle, L. Del Debbio, E. Kerrane, and J. Zanotti, Lattice Determination of the Hadronic 
Contribution to the Muon g — 2 using Dynamical Domain Wall Fermions, arXiv: 1107. 1497. 

G. Colangelo, S. Durr, A. Juttner, L. Lellouch, H. Leutwyler, et. al. Review of lattice results 
concerning low energy particle physics, Eur.Phys.J. C71 (2011) 1695, [arXiv: 1011.4408]. 

T. Blum, T. Doi, M. Hayakawa, T. Izubuchi, and N. Yamada, Determination of light quark 
masses from the electromagnetic splitting of pseudoscalar meson masses computed with two 
flavors of domain wall fermions, Phys. Rev. D76 (2007) 114508, [arXiv: 0708. 0484]. 

T. Blum et. at. Electromagnetic mass splittings of the low lying hadrons and quark masses from 
2-hl flavor lattice QCD-hQED, Phys. Rev. D82 (2010) 094508, [arXiv: 1006. 1311]. 



16 



[14] G. de Divitiis, P. Dimopoulos, R. Frezzotti, V. Lubicz, G. Martinelli, et. al., Isospin breaking 
effects due to the up-down mass difference in Lattice QCD, arXiv: 1110.6294. 

[15] A. Jiittner and M. Delia Morte, New ideas for g-2 on the lattice, PoS LAT2009 (2009) 143, 
[arXiv: 0910. 3755]. 

[16] M. Delia Morte and A. Jiittner, Quark disconnected diagrams in chiral perturbation theory, 
JHEPll (2010) 154, [arXiv: 1009. 3783]. 

[17] G. M. de Divitiis, R. Petronzio, and N. Tantalo, On the discretization of physical momenta in 
lattice QCD, Phys. Lett. B595 (2004) 408-413, [hep-lat/0405002]. 

[18] C. T. Sachrajda and G. Villadoro, Twisted boundary conditions in lattice simulations, Phys. 
Lett. B609 (2005) 73-85, [hep-lat/0411033]. 

[19] P. F. Bedaque and J.-W. Chen, Twisted valence quarks and hadron interactions on the lattice, 
Phys. Lett. B616 (2005) 208-214, [hep-lat/0412023]. 

[20] M. Delia Morte, B. Jager, A. Jiittner, and H. Wittig, The leading hadronic vacuum polarisation 
on the lattice, AIP Conf. Proc. 1343 (2011) 337-339, [arXiv: 1011 .5793]. 

[21] M. Delia Morte, B. Jager, A. Jiittner, and H. Wittig, Lattice Determination of the Anomalous 
Magnetic Moment of the Muon, PCS' Lattice 2011 (2011) 161, [arXiv: 1111 . 2193]. 

[22] E. de Rafael, Hadronic contributions to the muon g-2 and low-energy QCD, Phys. Lett. B322 
(1994) 239-246, [hep-pli/9311316]. 

[23] https://twiki.cern.ch/twiki/bin/view/CLS/WebIntro. 

[24] ALPHA Collaboration, K. Jansen and R. Sommer, 0(a) improvement of lattice QCD with two 
flavors of Wilson quarks, Nucl. Phys. B530 (1998) 185-203, [hep-lat/9803017]. 

[25] S. Capitani, M. Delia Morte, G. von Hippel, B. Knippschild, and H. Wittig, Scale setting via 
the Omega baryon mass, arXiv : 1110 . 6365. 

[26] L. Del Debbio, L. Giusti, M. Liischer, R. Petronzio, and N. Tantalo, QCD with light Wilson 
quarks on fine lattices (I): first experiences and physics results, JHEP 02 (2007) 056, 
[hep-lat/0610059]. 

[27] S. Capitani et. al., Mesonic and baryonic correlation functions at fine lattice spacings, PoS 
LAT2009 (2009) 095, [arXiv: 0910. 5578]. 

[28] B. B. Brandt et. al., Wilson fermions at fine lattice spacings: scale setting, pion form factors 
and (ff-2)^, PoS* LATTICE2010 (2010) 164, [arXiv: 1010.2390]. 

[29] G. Von Hippel and B. Knippschild. Private communication. 

[30] C. Michael, Fitting correlated data, Phys. Rev. D49 (1994) 2616-2619, [hep-lat/9310026]. 

[31] M. Liischer, Volume Dependence of the Energy Spectrum in Massive Quantum Field Theories. 

1. Stable Particle States, Commun.Math.Phys. 104 (1986) 177. 

[32] M. Liischer, Volume Dependence of the Energy Spectrum in Massive Quantum Field Theories. 

2. Scattering States, Commun.Math.Phys. 105 (1986) 153-188. 

[33] K. G. Chetyrkin, J. H. Kiihn, and M. Steinhauser, Three-loop polarization function and O(oij.) 
corrections to the production of heavy quarks, Nucl. Phys. B482 (1996) 213-240, 
[hep-ph/9606230]. 



17 



[34] ALPHA Collaboration, M. Delia Morte et. al, Computation of the strong coupling in QCD 
with two dynamical flavours, Nucl. Phys. B713 (2005) 378-406, [liep-lat/0411025]. 

[35] ALPHA Collaboration, M. Delia Morte et. al., Non-perturbative quark mass renormalization 
in two-flavor QCD, Nucl. Phys. B729 (2005) 117-134, [hep-lat/0507035]. 

[36] M. Delia Morte, R. Sommer, and S. Takeda, On cutoff effects in lattice QCD from short to long 
distances, Phys. Lett. B672 (2009) 407-412, [arXiv: 0807. 1120]. 

[37] P. Fritzsch, J. Heitger, and N. Tantalo, Non-perturbative improvement of quark mass 
renormalization in two-flavour lattice QCD, J HEP OH (2010) 074, [arXiv : 1004 . 3978]. 

[38] K. Nakamura and others (Particle Data Group), The review of particle physics, J. Phys. G 37 
(2010) 075021. 

[39] ETM Collaboration, D. B. Renner, X. Feng, K. Jansen, and M. Petschlies, Leading order 
hadronic contribution to g-2 from twisted mass QCD, PoS LATTICE2010 (2010) 155, 
[arXiv: 1011.4231]. 

[40] D. Bernecker and H. B. Meyer, Vector Correlators in Lattice QCD: methods and applications, 
arXiv: 1107.4388. 

[41] P. Boyle, J. Flynn, A. Jiittner, C. Sachrajda, and J. Zanotti, Hadronic form factors in Lattice 
QCD at small and vanishing momentum transfer, JHEP 0705 (2007) 016, [hep-lat/0703005]. 

[42] C. Thron, S. Dong, K. Liu, and H. Ying, Pade - Z(2) estimator of determinants, Phys. Rev. 
D57 (1998) 1642, [hep-lat/9707001]. 

[43] J. Foley, K. Jimmy Juge, A. O'Cais, M. Peardon, S. M. Ryan, et. al.. Practical all-to-all 

propagators for lattice QCD, Comput.Phys.Commun. 172 (2005) 145-162, [liep-lat/0505023]. 

[44] S. Collins, G. Bali, and A. Schafer, Disconnected contributions to hadronic structure: a new 
method for stochastic noise reduction, PoS" LAT2007 (2007) 141, [arXiv: 0709. 3217]. 

[45] C. Ehmann and G. S. Bali, Investigation of the rj' — rjc - mixing with improved stochastic 
estimators, Po^ LATTICE2008 (2008) 114, [arXiv : 0903 . 2947]. 

[46] V. Giilpers, Improved stochastic estimators for disconnected contributions to mesonic form 
factors in Lattice QCD, Diploma Thesis. 

littp://www.kph.uni-mainz.de/T/pub/diplonia/DipLTh_G uelpers.pdf. 



18 



